#define CATCH_CONFIG_MAIN

#include "catch.hpp"
#include "ENDFtk/section/2/151.hpp"

// other includes

// convenience typedefs
using namespace njoy::ENDFtk;
using Isotope = section::Type< 2, 151 >::Isotope;
using SpecialCase = section::Type< 2, 151 >::SpecialCase;
using SingleLevelBreitWigner = section::Type< 2, 151 >::SingleLevelBreitWigner;
using MultiLevelBreitWigner = section::Type< 2, 151 >::MultiLevelBreitWigner;
using ReichMoore = section::Type< 2, 151 >::ReichMoore;
using RMatrixLimited = section::Type< 2, 151 >::RMatrixLimited;
using UnresolvedEnergyIndependent = section::Type< 2, 151 >::UnresolvedEnergyIndependent;
using UnresolvedEnergyDependent = section::Type< 2, 151 >::UnresolvedEnergyDependent;
using UnresolvedEnergyDependentFissionWidths = section::Type< 2, 151 >::UnresolvedEnergyDependentFissionWidths;

SCENARIO( "Isotope" ){
  WHEN( "Negative NER" ){
    std::string ENDF =
      " 1.001000+3 9.991673-1          0          0         -1          0 125 2151    1\n";

    auto begin = ENDF.begin();
    auto end = ENDF.end();
    long lineNumber = 0;

    int MAT = 125;
    int MF = 2;
    int MT = 151;

    CHECK_THROWS(
      Isotope( begin, end, lineNumber, MAT, MF, MT )
    );
  }

  WHEN( "LRU=0" ){
    std::string ENDF =
      " 1.001000+3 1.000000+0          0          0          1          0 125 2151    2\n"
      " 1.000000-5 1.000000+5          0          0          0          0 125 2151    3\n"
      " 5.000000-1 1.276553+0          0          0          0          0 125 2151    4\n";

    auto begin = ENDF.begin();
    auto end = ENDF.end();
    long lineNumber = 0;

    int MAT = 125;
    int MF = 2;
    int MT = 151;

    THEN( "a SpecialCase can be constructed" ){
      Isotope iso( begin, end, lineNumber, MAT, MF, MT );

      CHECK( 1 == iso.resonanceRanges().size() );

      const auto& range = iso.resonanceRanges().back();
      CHECK( 1E-5 == Approx( range.EL() ) );
      CHECK( 1E-5 == Approx( range.lowerEnergy() ) );
      CHECK( 1E+5 == Approx( range.EH() ) );
      CHECK( 1E+5 == Approx( range.upperEnergy() ) );

      CHECK( 0 == range.LRU() );
      CHECK( 0 == range.type() );
      CHECK( 0 == range.LRF() );
      CHECK( 0 == range.representation() );
      CHECK( 0 == range.NRO() );
      CHECK( 0 == range.energyDependentScatteringRadius() );
      CHECK( 0 == range.NAPS() );
      CHECK( 0 == range.scatteringRadiusCalculationOption() );

      CHECK_NOTHROW( std::get< SpecialCase >
                       ( range.parameters() ) );
    }
  }
  WHEN( "LRU=1" ){
    WHEN( "there is only one SLBW range in the Isotope" ){
      std::string ENDF =
        " 6.114800E4 1.000000+0          0          0          1          06153 2151    2\n"
        " 1.000000-5 1.000000+0          1          1          0          06153 2151    3\n"
        " 6.000000+0 6.380900-1          0          0          1          06153 2151    4\n"
        " 1.466470+2 0.000000+0          0          0          6          16153 2151    5\n"
        " 1.690000-1 6.000000+0 8.018800-2 5.880000-4 7.960000-2 0.000000+06153 2151    6\n";

      auto begin = ENDF.begin();
      auto end = ENDF.end();
      long lineNumber = 0;
      int MAT = 6153;
      int MF = 2;
      int MT = 151;

      Isotope iso( begin, end, lineNumber, MAT, MF, MT );

      CHECK( 1 == iso.resonanceRanges().size() );

      const auto& range = iso.resonanceRanges().back();
      CHECK( 1E-5 == Approx( range.EL() ) );
      CHECK( 1E-5 == Approx( range.lowerEnergy() ) );
      CHECK( 1. == Approx( range.EH() ) );
      CHECK( 1. == Approx( range.upperEnergy() ) );

      CHECK( 1 == range.LRU() );
      CHECK( 1 == range.type() );
      CHECK( 1 == range.LRF() );
      CHECK( 1 == range.representation() );
      CHECK( 0 == range.NRO() );
      CHECK( 0 == range.energyDependentScatteringRadius() );
      CHECK( 0 == range.NAPS() );
      CHECK( 0 == range.scatteringRadiusCalculationOption() );

      CHECK_NOTHROW( std::get< SingleLevelBreitWigner >
                       ( range.parameters() ) );
    }

    WHEN( "there is only one MLBW range in the Isotope" ){
      std::string ENDF =
        " 6.114800E4 1.000000+0          0          0          1          06153 2151    2\n"
        " 1.000000-5 1.000000+0          1          2          0          06153 2151    3\n"
        " 6.000000+0 6.380900-1          0          0          1          06153 2151    4\n"
        " 1.466470+2 0.000000+0          0          0          6          16153 2151    5\n"
        " 1.690000-1 6.000000+0 8.018800-2 5.880000-4 7.960000-2 0.000000+06153 2151    6\n";

      auto begin = ENDF.begin();
      auto end = ENDF.end();
      long lineNumber = 0;
      int MAT = 6153;
      int MF = 2;
      int MT = 151;

      Isotope iso( begin, end, lineNumber, MAT, MF, MT );

      CHECK( 1 == iso.resonanceRanges().size() );

      const auto& range = iso.resonanceRanges().back();
      CHECK( 1E-5 == Approx( range.EL() ) );
      CHECK( 1E-5 == Approx( range.lowerEnergy() ) );
      CHECK( 1. == Approx( range.EH() ) );
      CHECK( 1. == Approx( range.upperEnergy() ) );

      CHECK( 1 == range.LRU() );
      CHECK( 1 == range.type() );
      CHECK( 2 == range.LRF() );
      CHECK( 2 == range.representation() );
      CHECK( 0 == range.NRO() );
      CHECK( 0 == range.energyDependentScatteringRadius() );
      CHECK( 0 == range.NAPS() );
      CHECK( 0 == range.scatteringRadiusCalculationOption() );

      CHECK_NOTHROW( std::get< MultiLevelBreitWigner >
                       ( range.parameters() ) );
    }

    WHEN( "there is only one Reich-Moore range in the Isotope" ){
      std::string ENDF =
        " 6.114800E4 1.000000+0          0          0          1          06153 2151    2\n"
        " 1.000000-5 1.000000+0          1          3          0          06153 2151    3\n"
        " 6.000000+0 6.380900-1          0          0          1          16153 2151    4\n"
        " 1.466470+2 0.000000+0          0          0          6          16153 2151    5\n"
        " 1.690000-1 6.000000+0 8.018800-2 5.880000-4 7.960000-2 0.000000+06153 2151    6\n";

      auto begin = ENDF.begin();
      auto end = ENDF.end();
      long lineNumber = 0;
      int MAT = 6153;
      int MF = 2;
      int MT = 151;

      Isotope iso( begin, end, lineNumber, MAT, MF, MT );

      CHECK( 1 == iso.resonanceRanges().size() );

      const auto& range = iso.resonanceRanges().back();
      CHECK( 1E-5 == Approx( range.EL() ) );
      CHECK( 1E-5 == Approx( range.lowerEnergy() ) );
      CHECK( 1. == Approx( range.EH() ) );
      CHECK( 1. == Approx( range.upperEnergy() ) );

      CHECK( 1 == range.LRU() );
      CHECK( 1 == range.type() );
      CHECK( 3 == range.LRF() );
      CHECK( 3 == range.representation() );
      CHECK( 0 == range.NRO() );
      CHECK( 0 == range.energyDependentScatteringRadius() );
      CHECK( 0 == range.NAPS() );
      CHECK( 0 == range.scatteringRadiusCalculationOption() );

      CHECK_NOTHROW( std::get< ReichMoore >
                       ( range.parameters() ) );
    }

    WHEN( "there is only one Adler-Adler range in the Isotope" ){
      std::string ENDF =
        " 6.114800E4 1.000000+0          0          0          1          06153 2151    2\n"
        " 1.000000-5 1.000000+0          1          4          0          06153 2151    3\n";

      auto begin = ENDF.begin();
      auto end = ENDF.end();
      long lineNumber = 0;
      int MAT = 6153;
      int MF = 2;
      int MT = 151;

      CHECK_THROWS(
        Isotope( begin, end, lineNumber, MAT, MF, MT )
      );
    }

    WHEN( "there is one Rmatrix-Limited range in the Isotope" ){
      std::string ENDF =
//    " 2.605400+4 5.347624+1          0          0          1          02625 2151     \n"  // HEAD
      " 2.605400+4 1.000000+0          0          0          1          02625 2151     \n"  // CONT (isotope)
      " 1.000000-5 1.036000+6          1          7          0          12625 2151     \n"  // CONT (range)
      " 0.000000+0 0.000000+0          0          3          5          02625 2151     \n"
      " 0.000000+0 0.000000+0          2          0         24          42625 2151     \n"
      " 0.000000+0 5.446635+1 0.000000+0 2.600000+1 1.000000+0 0.000000+02625 2151     \n"
      " 0.000000+0 0.000000+0 0.000000+0 1.020000+2 0.000000+0 0.000000+02625 2151     \n"
      " 1.000000+0 5.347624+1 0.000000+0 2.600000+1 5.000000-1 0.000000+02625 2151     \n"
      " 0.000000+0 1.000000+0 0.000000+0 2.000000+0 0.000000+0 1.000000+02625 2151     \n"
      " 5.000000-1 0.000000+0          0          0         12          22625 2151     \n"
      " 1.000000+0 0.000000+0 0.000000+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.000000+0 0.000000+0 5.000000-1 0.000000+0 5.437300-1 5.437300-12625 2151     \n"
      " 0.000000+0 0.000000+0          0        148        888        1482625 2151     \n"
      "-1.223300+6 1.000000+0 9.611086+5 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      "-2.500000+4 1.798017+0 1.032662+4 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.788000+3 1.455000+0 1.187354+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.287200+4 2.000000+0 2.000345+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.190500+4 2.000000+0 1.781791+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.881300+4 2.000000+0 5.551680+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.301700+5 3.000000+0 3.109384+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.475400+5 2.600000+0 3.083417+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.743000+5 2.000000+0 4.114302+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.887000+5 2.000000+0 3.684975+4 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.067800+5 4.000000-1 1.329600+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.236200+5 6.000000-1 8.984901+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.441500+5 2.000000+0 1.990200+4 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.573900+5 4.000000-1 4.010800+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.922000+5 1.100000+0 9.317101+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.072400+5 4.000000-1 6.319500+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.201400+5 4.000000-1 2.220800+4 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.316300+5 4.000000-1 3.519400+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.637100+5 4.000000-1 3.811500+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.712300+5 4.000000-1 4.752100+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.100200+5 4.000000-1 1.010700+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.112100+5 4.000000-1 3.435300+4 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.254800+5 4.000000-1 1.463000+4 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.302300+5 4.000000-1 3.839600+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.556500+5 4.000000-1 6.794200+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.832300+5 4.000000-1 2.227600+4 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.852200+5 4.000000-1 4.562300+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.041979+5 4.000000-1 3.716947+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.106047+5 4.000000-1 2.039219+4 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.268393+5 4.000000-1 9.906387+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.413763+5 4.000000-1 7.350561+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.587263+5 4.000000-1 6.587169+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.606915+5 4.000000-1 3.444498+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.928310+5 4.000000-1 2.908803+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.089289+5 4.000000-1 2.814805+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.420913+5 4.000000-1 1.435638+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.612270+5 4.000000-1 7.269890+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.726386+5 4.000000-1 2.202403+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.893352+5 4.000000-1 7.337199+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.056951+5 4.000000-1 1.145771-2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.197259+5 4.000000-1 5.241748+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.365430+5 4.000000-1 1.096020+4 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.484317+5 4.000000-1 1.583502+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.625934+5 4.000000-1 2.142637+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.934283+5 4.000000-1 3.224649+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.155311+5 4.000000-1 2.624618+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.191293+5 4.000000-1 3.436016+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.338111+5 4.000000-1 1.988010+4 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.440180+5 4.000000-1 7.598592+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.610880+5 4.000000-1 3.754241+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.712498+5 4.000000-1 4.100152+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.757820+5 4.000000-1 1.118407+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.770124+5 4.000000-1 1.126387-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.945130+5 4.000000-1 1.872254+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.005395+5 4.000000-1 9.809761-2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.045392+5 4.000000-1 4.966268+4 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.208361+5 4.000000-1 1.102448+4 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.270522+5 4.000000-1 5.401460+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.521109+5 4.000000-1 1.430473+4 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.601766+5 4.000000-1 3.057673+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.746035+5 4.000000-1 1.702876+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.836143+5 4.000000-1 3.736347-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.925091+5 4.000000-1 2.168022+4 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.989940+5 4.000000-1 3.905539+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.006835+6 7.000000-1 3.493911+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.007882+6 7.000000-1 1.857533+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.007962+6 7.000000-1 1.261542+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.012074+6 7.000000-1 4.408087+4 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.015771+6 7.000000-1 1.211998+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.016541+6 7.000000-1 7.275909+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.016807+6 7.000000-1 8.587253+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.023773+6 4.000000-1 1.382577+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.029100+6 4.000000-1 1.760486+4 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.030373+6 7.000000-1 1.878123+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.036413+6 4.000000-1 1.705277+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.036500+6 4.000000-1 1.196904+5 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.044895+6 7.000000-1 2.387808+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.048513+6 7.000000-1 2.865925+4 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.050702+6 7.000000-1 4.321734+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.053288+6 7.000000-1 1.668717+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.072431+6 7.000000-1 1.037008+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.073040+6 7.000000-1 4.756084+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.074304+6 7.000000-1 2.581614+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.076827+6 7.000000-1 2.561767+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.092165+6 7.000000-1 5.178618+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.094301+6 7.000000-1 4.002575+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.100703+6 7.000000-1 8.000866+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.101455+6 7.000000-1 1.629439+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.102090+6 7.000000-1 7.044857+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.105549+6 7.000000-1 2.826743+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.106805+6 7.000000-1 2.914327+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.113871+6 7.000000-1 6.341946+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.115013+6 7.000000-1 2.983193+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.116032+6 7.000000-1 2.011558+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.119116+6 7.000000-1 1.832122+4 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.125409+6 7.000000-1 1.602570+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.130170+6 7.000000-1 1.063606+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.137036+6 7.000000-1 1.177170+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.139917+6 7.000000-1 7.116193+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.148549+6 7.000000-1 6.479679+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.149079+6 7.000000-1 1.791819+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.153116+6 7.000000-1 1.291133+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.160226+6 7.000000-1 2.301283+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.166220+6 7.000000-1 1.626226+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.176926+6 7.000000-1 2.721416+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.177563+6 7.000000-1 2.833032+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.180217+6 7.000000-1 2.960697+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.184935+6 7.000000-1 1.349241+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.193010+6 7.000000-1 1.346723+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.211500+6 7.000000-1 8.129232+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.211900+6 7.000000-1 2.232713+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.216800+6 7.000000-1 7.532561+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.217600+6 7.000000-1 1.499205+4 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.219900+6 7.000000-1 1.368553+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.222600+6 7.000000-1 2.580568+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.228000+6 7.000000-1 1.045278+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.231200+6 7.000000-1 8.469479+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.235500+6 7.000000-1 1.296387+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.238200+6 7.000000-1 3.595342+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.261800+6 7.000000-1 5.152908+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.265000+6 7.000000-1 2.320658+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.270700+6 7.000000-1 5.740736+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.272800+6 7.000000-1 1.121145+4 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.277700+6 7.000000-1 1.963529+4 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.286400+6 7.000000-1 5.449450+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.288900+6 7.000000-1 4.610701+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.289300+6 7.000000-1 3.541559+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.294000+6 7.000000-1 1.513853+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.301464+6 7.000000-1 1.829046+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.303552+6 7.000000-1 6.868546+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.352640+6 7.000000-1 2.170071+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.382017+6 7.000000-1 8.166709-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.401250+6 7.000000-1 2.397060+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.407300+6 7.000000-1 1.597179+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.411100+6 7.000000-1 4.537008+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.412100+6 7.000000-1 4.967713+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.414000+6 7.000000-1 6.834148+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.434900+6 7.000000-1 1.972391+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.436800+6 7.000000-1 1.253945+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.459200+6 7.000000-1 4.672500+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.461300+6 7.000000-1 2.205592+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.472900+6 7.000000-1 4.017182+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.473500+6 7.000000-1 4.226838+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.475000+6 7.000000-1 5.630068+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.478600+6 7.000000-1 1.844172+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.501105+6 6.490600-1 1.598740+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.503308+6 6.490600-1 1.599913+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.505510+6 6.490600-1 1.601084+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      "-5.000000-1 0.000000+0          0          0         12          22625 2151     \n"
      " 1.000000+0 0.000000+0 0.000000+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.000000+0 1.000000+0 5.000000-1 0.000000+0 5.437300-1 5.437300-12625 2151     \n"
      " 0.000000+0 0.000000+0          0        176       1056        1762625 2151     \n"
      " 5.152000+4 3.600000-1 1.600200+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.359000+4 1.500000+0 1.700000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.545900+4 5.600000-1 3.200000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.432000+4 4.000000-1 2.601400+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.348000+4 5.900000-1 8.600000-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.135500+4 4.000000-1 6.601200+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.041600+5 4.000000-1 3.906100+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.158500+5 4.000000-1 2.152200+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.198300+5 4.000000-1 2.398400+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.263300+5 4.000000-1 6.520300+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.357600+5 3.500000-1 8.148600+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.504200+5 4.000000-1 7.445300+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.527000+5 4.000000-1 1.925500+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.571200+5 4.000000-1 1.077800+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.592800+5 4.000000-1 1.176400+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.652900+5 4.000000-1 6.101800+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.740300+5 4.000000-1 1.974800+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.779100+5 4.000000-1 7.004800+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.075000+5 4.000000-1 6.698100+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.076100+5 4.000000-1 1.182800+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.138500+5 4.000000-1 1.518300+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.253300+5 4.000000-1 4.758500+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.312300+5 5.900000-1 2.216667+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.377500+5 4.000000-1 7.063000+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.420400+5 4.000000-1 1.341400+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.623500+5 4.000000-1 4.649100+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.636000+5 4.000000-1 8.453800+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.714900+5 4.000000-1 7.560900+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.809900+5 3.200000-1 1.976800+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.893500+5 4.000000-1 6.920200+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.928600+5 4.000000-1 2.074600+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.035800+5 4.000000-1 1.919700+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.038200+5 4.000000-1 1.009800+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.091600+5 4.000000-1 2.961000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.117800+5 4.000000-1 4.457000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.260900+5 4.000000-1 9.192000+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.327500+5 4.000000-1 3.685100+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.367500+5 4.000000-1 5.497300+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.438500+5 4.000000-1 6.212900+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.573800+5 4.000000-1 2.848100+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.610300+5 5.900000-1 5.730000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.627300+5 5.900000-1 6.000000+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.673800+5 4.000000-1 2.910200+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.685200+5 4.000000-1 1.159600+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.701800+5 4.000000-1 1.031900+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.751300+5 4.000000-1 2.678900+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.847700+5 5.900000-1 7.800000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.896800+5 5.900000-1 1.240000+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.978500+5 5.900000-1 5.080000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.063100+5 4.000000-1 4.208800+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.163000+5 4.000000-1 7.167600+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.276600+5 4.000000-1 2.120400+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.277000+5 4.000000-1 3.344400+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.435300+5 4.000000-1 1.633100+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.464600+5 4.000000-1 8.464100+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.615900+5 4.000000-1 6.324700+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.686700+5 4.000000-1 5.382800+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.740000+5 4.000000-1 6.103500+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.906000+5 4.000000-1 1.008000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.930000+5 4.000000-1 1.442400+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.988200+5 4.000000-1 2.717300+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.109800+5 4.000000-1 1.087134+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.118797+5 4.000000-1 1.079896+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.242641+5 4.000000-1 4.204802+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.244715+5 4.000000-1 1.942133+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.260984+5 4.000000-1 1.099708+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.392393+5 4.000000-1 1.999082+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.490489+5 4.000000-1 6.418378+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.510297+5 4.000000-1 1.267160+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.612238+5 4.000000-1 9.983347+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.893758+5 4.000000-1 3.267754+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.064872+5 4.000000-1 3.363521+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.304033+5 4.000000-1 1.303078+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.361519+5 4.000000-1 8.801258+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.411656+5 4.000000-1 3.879805+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.569479+5 4.000000-1 4.685032+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.654921+5 4.000000-1 3.842323+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.848858+5 4.000000-1 1.630389-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.855506+5 4.000000-1 2.036212+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.079449+5 4.000000-1 3.262791+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.243508+5 4.000000-1 4.622010+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.390339+5 4.000000-1 5.374703+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.519040+5 4.000000-1 2.715538+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.913333+5 4.000000-1 9.720917+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.021403+5 4.000000-1 2.317511+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.052483+5 4.000000-1 2.727453+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.209221+5 4.000000-1 4.991260+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.230299+5 4.000000-1 1.540980-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.281600+5 4.000000-1 4.679765+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.388523+5 4.000000-1 6.492580+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.544406+5 4.000000-1 2.247298+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.742371+5 4.000000-1 1.183769+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.791545+5 4.000000-1 1.729120+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.851337+5 4.000000-1 2.502694+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.038761+5 4.000000-1 1.278049+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.148105+5 4.000000-1 5.179936+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.261085+5 4.000000-1 1.082181+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.458108+5 4.000000-1 1.290737+4 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.501513+5 4.000000-1 2.264682+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.561678+5 4.000000-1 2.404672+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.571613+5 4.000000-1 1.372520+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.706518+5 4.000000-1 3.140170+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.827969+5 4.000000-1 3.958198+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.918115+5 4.000000-1 4.349787-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.924344+5 4.000000-1 8.988154+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.943551+5 4.000000-1 5.400306+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.007936+6 7.000000-1 4.757654+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.010554+6 7.000000-1 1.788109+4 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.013662+6 7.000000-1 4.526361+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.014236+6 7.000000-1 1.639979+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.018384+6 7.000000-1 1.445604+4 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.020035+6 7.000000-1 2.270929+4 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.022489+6 7.000000-1 1.585767+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.024088+6 7.000000-1 2.741289+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.028337+6 7.000000-1 2.688911+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.028730+6 7.000000-1 1.169575+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.035029+6 7.000000-1 5.049021+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.041483+6 7.000000-1 6.932586+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.049194+6 7.000000-1 5.853234-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.051941+6 7.000000-1 2.732861-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.054551+6 7.000000-1 2.067607+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.058661+6 7.000000-1 7.427173+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.067529+6 7.000000-1 2.610806+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.079005+6 7.000000-1 4.984630+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.080066+6 7.000000-1 4.202750+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.087100+6 7.000000-1 2.263758+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.093357+6 7.000000-1 1.181581+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.095347+6 7.000000-1 3.341273+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.114263+6 7.000000-1 3.951291-3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.129215+6 7.000000-1 2.521005+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.136345+6 7.000000-1 3.742127+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.142586+6 7.000000-1 3.697277+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.142930+6 7.000000-1 9.441586+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.145446+6 7.000000-1 1.890100+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.153957+6 7.000000-1 6.507070+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.165404+6 7.000000-1 2.397530+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.169624+6 7.000000-1 2.521247+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.172316+6 7.000000-1 4.367571+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.178299+6 7.000000-1 2.281827+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.186150+6 7.000000-1 2.175313+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.196735+6 7.000000-1 1.551876+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.197633+6 7.000000-1 1.472436+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.200900+6 7.000000-1 1.622288+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.203200+6 7.000000-1 2.045252+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.203900+6 7.000000-1 2.976736+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.211100+6 7.000000-1 8.314176+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.214200+6 7.000000-1 7.927457+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.228800+6 7.000000-1 1.436400+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.235900+6 7.000000-1 3.275355+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.246100+6 7.000000-1 3.361403+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.250500+6 7.000000-1 2.404723+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.258500+6 7.000000-1 8.707712+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.269500+6 7.000000-1 3.799778+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.288500+6 7.000000-1 4.659112+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.291400+6 7.000000-1 1.021775+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.293500+6 7.000000-1 7.503388+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.294400+6 7.000000-1 1.665153+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.295600+6 7.000000-1 4.178903+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.299000+6 7.000000-1 4.747361+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.299900+6 7.000000-1 3.440912+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.321001+6 7.000000-1 7.322842+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.325927+6 7.000000-1 1.240116+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.341047+6 7.000000-1 7.350127+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.365848+6 7.000000-1 5.123859-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.376402+6 7.000000-1 1.255995+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.384840+6 7.000000-1 3.267945+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.389363+6 7.000000-1 2.240603+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.390975+6 7.000000-1 4.612233+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.394987+6 7.000000-1 5.526015+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.415000+6 7.000000-1 3.262509+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.433400+6 7.000000-1 8.731223-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.437300+6 7.000000-1 2.112829+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.440300+6 7.000000-1 1.097058+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.458700+6 7.000000-1 1.608010+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.462300+6 7.000000-1 2.553174+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.470900+6 7.000000-1 1.946542+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      "-1.500000+0 0.000000+0          0          0         12          22625 2151     \n"
      " 1.000000+0 0.000000+0 0.000000+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.000000+0 1.000000+0 5.000000-1 0.000000+0 5.437300-1 5.437300-12625 2151     \n"
      " 0.000000+0 0.000000+0          0        196       1176        1962625 2151     \n"
      " 3.099000+3 5.900000-1 1.400000-3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.358100+4 5.900000-1 1.750000-2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.927800+4 5.900000-1 2.750000-2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.303000+4 5.900000-1 2.675000-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.822000+4 5.900000-1 9.250000-2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.526000+4 5.900000-1 1.400000-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.844000+4 5.900000-1 1.400000+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.120000+4 5.900000-1 1.200000-2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.014000+4 5.900000-1 2.900000-2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.920000+4 5.900000-1 2.700000-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.876000+4 5.900000-1 1.500000-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.580000+4 4.000000-1 9.166800-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.720000+4 4.000000-1 2.743400+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.128000+4 5.900000-1 1.900000-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.315500+4 4.000000-1 3.876100+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.733000+4 5.900000-1 2.900000-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.135500+4 5.900000-1 1.300000-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.777000+4 5.900000-1 6.100000-2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.984000+4 5.900000-1 9.600000-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.017100+5 5.900000-1 1.900000-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.125800+5 5.900000-1 7.000000-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.129600+5 5.900000-1 3.300000-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.158500+5 5.500000-1 1.350000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.197800+5 5.300000-1 1.300000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.265300+5 8.000000-1 2.112100+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.379300+5 5.900000-1 6.300000-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.410100+5 5.900000-1 5.900000-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.455100+5 5.900000-1 4.100000-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.531900+5 4.000000-1 2.270600+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.572000+5 5.600000-1 1.085000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.645100+5 1.400000+0 5.649700+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.890900+5 5.900000-1 2.500000-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.944800+5 4.000000-1 4.153200+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.039700+5 4.000000-1 1.609000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.076000+5 5.900000-1 1.000000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.100100+5 5.900000-1 3.870000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.257100+5 5.900000-1 5.700000+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.312200+5 4.000000-1 1.590600+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.337100+5 5.900000-1 1.295000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.420600+5 5.900000-1 8.800000+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.452900+5 4.000000-1 2.665600+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.548900+5 5.900000-1 3.775000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.630200+5 4.000000-1 1.642900+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.829400+5 4.000000-1 8.750400+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.091400+5 5.900000-1 1.310000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.220200+5 4.000000-1 1.284700+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.326900+5 5.900000-1 1.455000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.452700+5 4.000000-1 1.177500+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.570900+5 4.000000-1 1.296200+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.610900+5 4.000000-1 2.322500+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.848400+5 4.000000-1 2.358900+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.897200+5 4.000000-1 5.620500+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.899100+5 4.000000-1 1.699400+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.979100+5 4.300000-1 2.239900+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.988000+5 4.000000-1 2.994400+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.123900+5 4.000000-1 9.687500+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.275300+5 5.900000-1 1.115000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.277000+5 5.900000-1 1.410000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.469700+5 4.000000-1 6.387000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.589700+5 5.900000-1 1.390000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.614600+5 5.900000-1 3.205000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.654500+5 5.900000-1 1.810000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.703300+5 4.000000-1 1.404700+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.786900+5 4.000000-1 8.111300+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.901500+5 4.000000-1 4.212400+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.980600+5 4.000000-1 1.854200+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.113799+5 4.000000-1 2.661896+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.242200+5 5.900000-1 7.000000+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.305283+5 4.000000-1 2.377581+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.525716+5 4.000000-1 1.806954+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.598908+5 4.000000-1 5.505216+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.651993+5 4.000000-1 7.960016+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.687365+5 4.000000-1 9.021971+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.809309+5 4.000000-1 4.985751+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.952300+5 5.900000-1 5.050000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.042492+5 4.000000-1 8.322616+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.044843+5 4.000000-1 4.213639+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.212139+5 4.000000-1 7.294993+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.385500+5 5.900000-1 2.755000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.470683+5 4.000000-1 5.606430+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.567200+5 5.900000-1 2.000000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.699057+5 4.000000-1 2.493408+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.768422+5 4.000000-1 3.262786+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.800487+5 4.000000-1 4.561021+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.805754+5 4.000000-1 3.858567+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.818110+5 4.000000-1 5.694376+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.835464+5 4.000000-1 4.575830+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.840877+5 4.000000-1 2.825082-4 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.876284+5 4.000000-1 4.352441+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.894095+5 4.000000-1 4.969923+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.978326+5 4.000000-1 6.647285+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.087790+5 4.000000-1 5.916424+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.226067+5 4.000000-1 1.205360+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.265519+5 4.000000-1 4.449364+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.387100+5 5.900000-1 2.675000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.494766+5 4.000000-1 8.185866+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.735575+5 4.000000-1 8.873871+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.909786+5 4.000000-1 6.627769+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.954090+5 4.000000-1 9.319217+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.070882+5 4.000000-1 3.330302+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.129077+5 4.000000-1 8.293181+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.186005+5 4.000000-1 9.917532+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.254062+5 4.000000-1 2.333022+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.290519+5 4.000000-1 6.492107+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.357682+5 4.000000-1 4.500824+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.412780+5 4.000000-1 1.254818+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.434657+5 4.000000-1 6.169918+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.474795+5 4.000000-1 2.071445+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.522371+5 4.000000-1 2.188346+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.591291+5 4.000000-1 4.349779+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.763763+5 4.000000-1 2.369700+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.881116+5 4.000000-1 1.688519+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.960101+5 4.000000-1 1.219338+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.977190+5 4.000000-1 1.028863+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.099892+5 4.000000-1 1.294740+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.144508+5 4.000000-1 2.595135+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.164981+5 4.000000-1 1.931471+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.187953+5 4.000000-1 1.313353+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.227527+5 4.000000-1 1.658686+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.273890+5 4.000000-1 2.096364+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.448935+5 4.000000-1 7.524011+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.636213+5 4.000000-1 1.307320+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.680020+5 4.000000-1 3.850489+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.784231+5 4.000000-1 8.367731+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.838438+5 4.000000-1 4.147539+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.900145+5 4.000000-1 5.543467-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.002368+6 7.000000-1 1.518759+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.005913+6 7.000000-1 7.659975+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.010327+6 7.000000-1 7.439483+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.012920+6 7.000000-1 7.671746+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.022815+6 7.000000-1 2.289961-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.026431+6 7.000000-1 1.268982+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.028429+6 7.000000-1 1.547053+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.031001+6 7.000000-1 3.223212+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.031465+6 7.000000-1 3.598004+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.032329+6 7.000000-1 4.183987+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.051001+6 7.000000-1 6.110995+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.056050+6 7.000000-1 4.184111+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.058038+6 7.000000-1 6.158646+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.062925+6 7.000000-1 2.761361+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.072620+6 7.000000-1 9.024670+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.083488+6 7.000000-1 1.088390+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.085559+6 7.000000-1 1.483430-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.091677+6 7.000000-1 1.194235+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.095462+6 7.000000-1 1.894608+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.099070+6 7.000000-1 5.528152+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.101233+6 7.000000-1 1.596818+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.104403+6 7.000000-1 1.760413+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.106613+6 7.000000-1 1.989903+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.107490+6 7.000000-1 1.274670+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.109121+6 7.000000-1 3.617540+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.110144+6 7.000000-1 2.370016+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.117903+6 7.000000-1 8.918474+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.118632+6 7.000000-1 7.202237+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.119608+6 7.000000-1 8.576966+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.125225+6 7.000000-1 3.283345+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.128490+6 7.000000-1 1.075546+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.135245+6 7.000000-1 6.675872-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.141235+6 7.000000-1 6.903190+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.144014+6 7.000000-1 8.474489+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.160850+6 7.000000-1 1.739419+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.166727+6 7.000000-1 3.483328+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.172094+6 7.000000-1 5.614840-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.179700+6 7.000000-1 1.049276+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.190647+6 7.000000-1 1.988563+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.191495+6 7.000000-1 9.753071-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.195039+6 7.000000-1 1.337952+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.197252+6 7.000000-1 2.097729+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.200500+6 7.000000-1 1.226610+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.214900+6 7.000000-1 6.750564+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.220700+6 7.000000-1 6.827222-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.274800+6 7.000000-1 5.396743+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.288100+6 7.000000-1 2.739994+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.297300+6 7.000000-1 1.083002+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.300387+6 7.000000-1 5.493265+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.306103+6 7.000000-1 1.778696+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.307309+6 7.000000-1 1.608504-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.312399+6 7.000000-1 6.039068+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.318355+6 7.000000-1 1.022869+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.324614+6 7.000000-1 7.662769+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.362008+6 7.000000-1 9.721002+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.370898+6 7.000000-1 1.373088+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.382499+6 7.000000-1 3.325686+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.396057+6 7.000000-1 1.046229+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.403500+6 7.000000-1 3.838244+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.407800+6 7.000000-1 8.793682+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.421700+6 7.000000-1 9.757445+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.430000+6 7.000000-1 9.268406-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.432400+6 7.000000-1 7.682903-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.439300+6 7.000000-1 1.981232+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.441300+6 7.000000-1 5.373680-3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.465800+6 7.000000-1 4.801503+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.467900+6 7.000000-1 8.062626+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.468900+6 7.000000-1 2.746135+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.475500+6 7.000000-1 2.636608+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.496700+6 7.000000-1 5.591693+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.500000+0 0.000000+0          0          0         12          22625 2151     \n"
      " 1.000000+0 0.000000+0 0.000000+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.000000+0 2.000000+0 5.000000-1 0.000000+0 5.437300-1 5.437300-12625 2151     \n"
      " 0.000000+0 0.000000+0          0        174       1044        1742625 2151     \n"
      " 9.480000+3 2.700000-1 1.200000+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.118000+4 3.500000-1 3.850100+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.445000+4 3.500000-1 7.000200-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.064000+4 4.800000-1 5.000100+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.914000+4 4.100000-1 8.500200-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.354000+4 3.000000-1 8.502400+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.535000+4 3.400000-1 1.600400+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.206900+5 4.000000-1 9.119800-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.207700+5 4.500000-1 1.613900+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.428200+5 9.600000-1 1.200000+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.646500+5 4.000000-1 2.848500+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.824200+5 5.010000-1 4.015000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.921200+5 9.600000-1 1.200000+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.977000+5 9.600000-1 7.600000-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.099900+5 4.000000-1 2.473800+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.257700+5 4.000000-1 1.283500+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.258000+5 4.000000-1 1.312500+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.280900+5 4.000000-1 3.677300+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.293000+5 4.000000-1 5.195600+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.336800+5 4.000000-1 1.611700+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.549000+5 4.000000-1 2.828300+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.629000+5 4.000000-1 1.811800+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.642500+5 4.000000-1 8.732600+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.773200+5 4.000000-1 1.202800+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.059800+5 4.000000-1 1.287800+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.241600+5 4.000000-1 6.035600+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.677300+5 4.000000-1 2.369800+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.830000+5 4.000000-1 2.006800+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.907600+5 4.000000-1 1.392500+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.911900+5 4.000000-1 3.733800+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 3.989500+5 4.000000-1 2.217800+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.029600+5 4.000000-1 1.621900+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.095000+5 4.000000-1 3.039900+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.136700+5 4.000000-1 1.566300+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.167000+5 4.000000-1 4.239100+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.279900+5 4.000000-1 4.708300+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.321100+5 4.000000-1 1.917200+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.428100+5 4.000000-1 1.435500+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.460300+5 4.000000-1 2.966900+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.465700+5 4.000000-1 5.843200+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.590400+5 4.000000-1 1.394400+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.655400+5 4.000000-1 2.570000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.760200+5 4.000000-1 8.034900+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.762900+5 4.000000-1 2.767100+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.820800+5 4.000000-1 1.243900+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 4.857700+5 4.000000-1 5.517600+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.073475+5 4.000000-1 1.379720+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.102716+5 4.000000-1 6.736567+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.104292+5 4.000000-1 1.103354+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.154764+5 4.000000-1 2.801142+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.190815+5 4.000000-1 1.666742+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.216600+5 4.000000-1 4.820773+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.341374+5 4.000000-1 5.966375+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.385437+5 4.000000-1 1.679649+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.496540+5 4.000000-1 1.547438+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.561857+5 4.000000-1 4.058745+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.694581+5 4.000000-1 1.073030+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.716764+5 4.000000-1 2.133730+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.761440+5 4.000000-1 2.209928+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.768834+5 4.000000-1 3.343362+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.791168+5 4.000000-1 5.318386+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.852562+5 4.000000-1 3.838181+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.856907+5 4.000000-1 4.050193+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.882995+5 4.000000-1 4.559138+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.953607+5 4.000000-1 4.823196+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.972367+5 4.000000-1 5.953381+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 5.990237+5 4.000000-1 3.782001+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.167820+5 4.000000-1 2.131430+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.186659+5 4.000000-1 4.343671+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.348638+5 4.000000-1 1.106123+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.389130+5 4.000000-1 1.426816+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.428880+5 4.000000-1 5.749459+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.435311+5 4.000000-1 1.272912-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.485029+5 4.000000-1 1.072359+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.529115+5 4.000000-1 5.154546+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.547732+5 4.000000-1 3.807383+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.631077+5 4.000000-1 2.946856+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.648455+5 4.000000-1 7.303448+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.749092+5 4.000000-1 5.093406+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.986631+5 4.000000-1 1.384974+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.048440+5 4.000000-1 1.236817+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.172003+5 4.000000-1 1.853490+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.310510+5 4.000000-1 9.744808+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.415731+5 4.000000-1 5.859396+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.446768+5 4.000000-1 1.510914+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.537304+5 4.000000-1 3.635781+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.606821+5 4.000000-1 7.738128+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.695500+5 9.600000-1 2.750000+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.755800+5 9.600000-1 2.330000+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.839034+5 4.000000-1 6.394942+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.883167+5 4.000000-1 9.858968-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.070300+5 9.600000-1 7.950000+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.199535+5 4.000000-1 5.422480+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.622708+5 4.000000-1 2.318639+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.723056+5 4.000000-1 7.052260-2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.238092+5 4.000000-1 2.697234-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.329524+5 4.000000-1 1.072767+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.365564+5 4.000000-1 1.358183+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.616729+5 4.000000-1 1.123392+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.911457+5 4.000000-1 9.450324+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.020395+6 7.000000-1 1.286129+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.025416+6 7.000000-1 2.147429+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.032867+6 7.000000-1 8.407657-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.034654+6 7.000000-1 1.485668+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.036976+6 7.000000-1 2.008344+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.037442+6 7.000000-1 1.620103+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.044659+6 7.000000-1 2.324646-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.056393+6 7.000000-1 1.582706+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.063226+6 7.000000-1 1.769481+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.068528+6 7.000000-1 2.381345+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.069529+6 7.000000-1 2.270916+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.071085+6 7.000000-1 5.606077+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.083799+6 7.000000-1 6.748964+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.085802+6 7.000000-1 1.336727+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.087257+6 7.000000-1 4.625587+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.089521+6 7.000000-1 4.606651-2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.090365+6 7.000000-1 2.133190+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.097836+6 7.000000-1 4.604364-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.107098+6 7.000000-1 1.165450+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.110852+6 7.000000-1 3.184407-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.138229+6 7.000000-1 1.098227+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.138864+6 7.000000-1 2.995481+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.140031+6 7.000000-1 1.256156+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.144947+6 7.000000-1 5.943933-3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.160565+6 7.000000-1 1.941732+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.164436+6 7.000000-1 4.292274+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.172409+6 7.000000-1 2.776853+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.182174+6 7.000000-1 2.143896+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.183615+6 7.000000-1 2.970873+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.188969+6 7.000000-1 2.718275+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.193005+6 7.000000-1 1.868313+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.200100+6 7.000000-1 1.243752+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.206900+6 7.000000-1 1.720366+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.218000+6 7.000000-1 1.068909+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.226500+6 7.000000-1 1.388803+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.226900+6 7.000000-1 6.263633+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.231900+6 7.000000-1 4.681152-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.234700+6 7.000000-1 2.314244+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.242600+6 7.000000-1 2.928929+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.243000+6 7.000000-1 4.564524+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.245700+6 7.000000-1 3.598557-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.246500+6 7.000000-1 6.446429+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.247300+6 7.000000-1 1.382621+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.249300+6 7.000000-1 8.079770+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.250100+6 7.000000-1 8.699172+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.252900+6 7.000000-1 1.340302+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.255300+6 7.000000-1 3.059559-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.255700+6 7.000000-1 3.350943+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.258100+6 7.000000-1 6.970356+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.262200+6 7.000000-1 3.769258+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.264600+6 7.000000-1 1.110896+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.271100+6 7.000000-1 3.214584+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.272400+6 7.000000-1 2.032512+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.280200+6 7.000000-1 6.998520+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.285600+6 7.000000-1 8.680330+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.291900+6 7.000000-1 2.057113+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.296100+6 7.000000-1 3.140329+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.310522+6 7.000000-1 1.626881+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.342597+6 7.000000-1 2.866731+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.353687+6 7.000000-1 6.316532+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.373551+6 7.000000-1 1.034589-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.386529+6 7.000000-1 2.898360+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.392287+6 7.000000-1 2.060996+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.400200+6 7.000000-1 1.115285+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.404400+6 7.000000-1 1.298045+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.419800+6 7.000000-1 1.121747+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.432900+6 7.000000-1 5.891910-2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.435400+6 7.000000-1 5.562327-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.438800+6 7.000000-1 2.492189+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.455700+6 7.000000-1 2.276831+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.459800+6 7.000000-1 2.752649+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.465300+6 7.000000-1 1.902441+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.469400+6 7.000000-1 2.018815+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.480100+6 7.000000-1 1.766491+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.500000+0 0.000000+0          0          0         12          22625 2151     \n"
      " 1.000000+0 0.000000+0 0.000000+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.000000+0 2.000000+0 5.000000-1 0.000000+0 5.437300-1 5.437300-12625 2151     \n"
      " 0.000000+0 0.000000+0          0         43        258         432625 2151     \n"
      " 1.264000+5 1.100000+0 2.900000+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.504700+5 9.600000-1 2.600000+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.779400+5 9.600000-1 1.400000+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.823600+5 4.000000-1 5.609500+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.945200+5 9.600000-1 5.300000-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 2.635000+5 4.000000-1 3.584800+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 6.241610+5 4.000000-1 1.280927+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.151255+5 4.000000-1 3.506706+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.196564+5 4.000000-1 4.798990+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.282568+5 4.000000-1 6.731346+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.406592+5 4.000000-1 1.917303+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.506336+5 4.000000-1 3.794667+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.555251+5 4.000000-1 2.078499+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.643228+5 4.000000-1 1.046145+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.671751+5 4.000000-1 1.961934+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.698952+5 4.000000-1 1.742457+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.759082+5 4.000000-1 1.815106+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.823973+5 4.000000-1 3.856682+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 7.965074+5 4.000000-1 6.345531+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.122098+5 4.000000-1 2.518495-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.307875+5 4.000000-1 7.280372+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.465205+5 4.000000-1 1.165599+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.564689+5 4.000000-1 1.551788+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.630488+5 4.000000-1 1.511368+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.699976+5 4.000000-1 1.737354+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.757905+5 4.000000-1 9.092705+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.814201+5 4.000000-1 1.593016+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 8.864225+5 4.000000-1 3.976289+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.012799+5 4.000000-1 1.208135+3 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.123647+5 4.000000-1 3.051593+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.382464+5 4.000000-1 5.450844+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.393989+5 4.000000-1 8.281391+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.413543+5 4.000000-1 1.625061+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.422685+5 4.000000-1 1.445374+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.430822+5 4.000000-1 1.347471+0 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.527008+5 4.000000-1 4.198989-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.555271+5 4.000000-1 1.052502+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.565799+5 4.000000-1 6.338464+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.584438+5 4.000000-1 1.034887+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.769656+5 4.000000-1 2.873313+2 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.794482+5 4.000000-1 8.738206+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 9.805248+5 4.000000-1 9.788140-1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 1.001152+6 4.000000-1 6.940802+1 0.000000+0 0.000000+0 0.000000+02625 2151     \n"
      " 0.000000+0 0.000000+0          0          0          0          02625 2  0     \n";

      auto begin = ENDF.begin();
      auto end = ENDF.end();
      long lineNumber = 0;
      int MAT = 2625;
      int MF = 2;
      int MT = 151;

      Isotope iso( begin, end, lineNumber, MAT, MF, MT );

      CHECK( 1 == iso.resonanceRanges().size() );

      const auto& range = iso.resonanceRanges().back();
      CHECK( 1E-5 == Approx( range.EL() ) );
      CHECK( 1E-5 == Approx( range.lowerEnergy() ) );
      CHECK( 1.036e+6 == Approx( range.EH() ) );
      CHECK( 1.036e+6 == Approx( range.upperEnergy() ) );

      CHECK( 1 == range.LRU() );
      CHECK( 1 == range.type() );
      CHECK( 7 == range.LRF() );
      CHECK( 7 == range.representation() );
      CHECK( 0 == range.NRO() );
      CHECK( 0 == range.energyDependentScatteringRadius() );
      CHECK( 1 == range.NAPS() );
      CHECK( 1 == range.scatteringRadiusCalculationOption() );

      CHECK_NOTHROW(
        std::get< RMatrixLimited >(
            range.parameters() )
      );
    }

    WHEN( "there is more than one range in the Isotope" ){
      std::string ENDF =
        " 6.114800E4 1.000000+0          0          0          2          06153 2151    2\n"
        " 1.000000-5 1.000000+0          1          1          0          06153 2151    3\n"
        " 6.000000+0 6.380900-1          0          0          1          06153 2151    4\n"
        " 1.466470+2 0.000000+0          0          0          6          16153 2151    5\n"
        " 1.690000-1 6.000000+0 8.018800-2 5.880000-4 7.960000-2 0.000000+06153 2151    6\n"
        " 1.000000-5 1.000000+0          1          3          0          06153 2151    7\n"
        " 6.000000+0 6.380900-1          0          0          1          16153 2151    8\n"
        " 1.466470+2 0.000000+0          0          0          6          16153 2151    9\n"
        " 1.690000-1 6.000000+0 8.018800-2 5.880000-4 7.960000-2 0.000000+06153 2151   10\n";

      auto begin = ENDF.begin();
      auto end = ENDF.end();
      long lineNumber = 0;
      int MAT = 6153;
      int MF = 2;
      int MT = 151;

      Isotope iso( begin, end, lineNumber, MAT, MF, MT );
      CHECK( 2 == iso.resonanceRanges().size() );

      CHECK_NOTHROW(
        std::get< SingleLevelBreitWigner >
                       ( iso.resonanceRanges().front().parameters() ) );
      CHECK_NOTHROW(
          std::get< ReichMoore >
                       ( iso.resonanceRanges().back().parameters() ) );
    }

    WHEN( "invalid LRF values" ){
      njoy::Log::info( "Testing invalid LRF with LRU=1" );
      int MAT = 125;
      int MF = 2;
      int MT = 151;
      THEN( "LRF=0, an exception is thrown" ){
        std::string ENDF =
          " 1.001000+3 1.000000+0          0          0          1          0 125 2151    2\n"
          " 1.000000-5 1.000000+5          1          0          0          0 125 2151    3\n";
        auto begin = ENDF.begin();
        auto end = ENDF.end();
        long lineNumber = 0;

        CHECK_THROWS(
          Isotope( begin, end, lineNumber, MAT, MF, MT ) );
      }
      THEN( "LRF=5, an exception is thrown" ){
        std::string ENDF =
          " 1.001000+3 1.000000+0          0          0          1          0 125 2151    2\n"
          " 1.000000-5 1.000000+5          1          5          0          0 125 2151    3\n";
        auto begin = ENDF.begin();
        auto end = ENDF.end();
        long lineNumber = 0;

        CHECK_THROWS(
          Isotope( begin, end, lineNumber, MAT, MF, MT ) );
      }
      THEN( "LRF=6, an exception is thrown" ){
        std::string ENDF =
          " 1.001000+3 1.000000+0          0          0          1          0 125 2151    2\n"
          " 1.000000-5 1.000000+5          1          6          0          0 125 2151    3\n";
        auto begin = ENDF.begin();
        auto end = ENDF.end();
        long lineNumber = 0;

        CHECK_THROWS(
          Isotope( begin, end, lineNumber, MAT, MF, MT ) );
      }

      THEN( "LRF=8, an exception is thrown" ){
        std::string ENDF =
          " 1.001000+3 1.000000+0          0          0          1          0 125 2151    2\n"
          " 1.000000-5 1.000000+5          1          8          0          0 125 2151    3\n";
        auto begin = ENDF.begin();
        auto end = ENDF.end();
        long lineNumber = 0;

        CHECK_THROWS(
          Isotope( begin, end, lineNumber, MAT, MF, MT ) );
      }
    }
  }
  WHEN( "LRU=2" ){
    WHEN( "LRF=1" ){
      WHEN( "LFW=0, energy-independent" ){
        std::string ENDF =
          " 6.816700+4 1.000000+0          0          0          1          05655 2151\n"
          // base
          " 2.300000+4 1.000000+5          2          1          0          05655 2151\n"
          // range CONT
          " 0.000000+0 6.233000-1          0          0          3          05655 2151\n"
          // L=0 LIST
          " 1.387090+2 0.000000+0          0          0          6          15655 2151\n"
          " 4.400000+3 5.000000-1 1.000000+0 4.400000-1 5.000000-2 0.000000+05655 2151\n"
          // L=1 LIST
          " 1.387090+2 0.000000+0          1          0         12          25655 2151\n"
          " 4.400000+3 5.000000-1 1.000000+0 5.280000-1 9.000000-2 0.000000+05655 2151\n"
          " 2.200000+3 1.500000+0 1.000000+0 2.640000-1 9.000000-2 0.000000+05655 2151\n"
          // L=2 LIST
          " 1.387090+2 0.000000+0          2          0         12          25655 2151\n"
          " 2.200000+3 1.500000+0 1.000000+0 3.300000-2 5.000000-2 0.000000+05655 2151\n"
          " 1.466670+3 2.500000+0 1.000000+0 2.200000-2 5.000000-2 0.000000+05655 2151\n"
          " 0.000000+0 0.000000+0          0          0          0          05655 2  0\n"
          " 0.000000+0 0.000000+0          0          0          0          05655 0  0\n";

        auto begin = ENDF.begin();
        auto end = ENDF.end();
        long lineNumber = 0;

        int MAT = 5655;
        int MF = 2;
        int MT = 151;

        THEN( "the parameters can be extracted and verified" ){
          Isotope iso( begin, end, lineNumber, MAT, MF, MT );
          CHECK_NOTHROW(
            std::get<
              UnresolvedEnergyIndependent >
              ( iso.resonanceRanges().front().parameters() )
          );
        }
      }
      WHEN( "LFW=1, fission widths given" ){
        std::string ENDF =
          " 6.816700+4 1.000000+0          0          1          1          09440 2151\n"
          " 5.700000+3 4.000000+4          2          1          0          09440 2151\n"
          // LIST
          " 0.000000+0 8.880000-1          0          0         14          39440 2151\n"
          " 5.700000+3 7.000000+3 8.000000+3 9.000000+3 1.000000+4 1.200000+49440 2151\n"
          " 1.400000+4 1.600000+4 1.800000+4 2.000000+4 2.500000+4 3.000000+49440 2151\n"
          " 3.500000+4 4.000000+4                                            9440 2151\n"
          // L=0 CONT
          " 2.379920+2 0.000000+0          0          0          1          09440 2151\n"
          // L=0 LIST
          " 0.000000+0 0.000000+0          0          1         20          09440 2151\n"
          " 1.310000+1 5.000000-1 1.000000+0 1.572000-3 3.100000-2 0.000000+09440 2151\n"
          " 1.256000-3 1.544000-3 1.825000-3 2.025000-3 2.119000-3 2.051000-39440 2151\n"
          " 1.992000-3 1.879000-3 1.860000-3 1.838000-3 1.694000-3 1.581000-39440 2151\n"
          " 1.481000-3 1.403000-3                                            9440 2151\n"
          // L=1 CONT
          " 2.379920+2 0.000000+0          1          0          2          09440 2151\n"
          // L=1,J=0 LIST
          " 0.000000+0 0.000000+0          1          1         20          09440 2151\n"
          " 1.310000+1 5.000000-1 1.000000+0 3.013000-3 3.100000-2 0.000000+09440 2151\n"
          " 4.314000-3 4.572000-3 4.740000-3 5.000000-3 5.520000-3 7.057000-39440 2151\n"
          " 8.251000-3 9.276000-3 9.930000-3 1.035000-2 1.210000-2 1.341000-29440 2151\n"
          " 1.456000-2 1.542000-2                                            9440 2151\n"
          // L=1,J=1 LIST
          " 0.000000+0 0.000000+0          1          2         20          09440 2151\n"
          " 6.697000+0 1.500000+0 1.000000+0 1.540000-3 3.100000-2 0.000000+09440 2151\n"
          " 1.207000-3 1.299000-3 1.326000-3 1.397000-3 1.544000-3 1.973000-39440 2151\n"
          " 2.317000-3 2.581000-3 2.811000-3 3.075000-3 3.494000-3 3.887000-39440 2151\n"
          " 4.290000-3 4.586000-3                                            9440 2151\n"
          // L=2 CONT
          " 2.379920+2 0.000000+0          2          0          2          09440 2151\n"
          // L=2,J=0 LIST
          " 0.000000+0 0.000000+0          2          1         20          09440 2151\n"
          " 6.697000+0 1.500000+0 1.000000+0 8.304300-4 3.100000-2 0.000000+09440 2151\n"
          " 1.000000-9 1.000000-9 1.000000-9 1.000000-9 1.000000-9 1.000000-99440 2151\n"
          " 1.000000-9 1.000000-9 1.000000-9 1.000000-9 1.000000-9 1.000000-99440 2151\n"
          " 1.000000-9 1.000000-9                                            9440 2151\n"
          // L=2,J=1 LIST
          " 0.000000+0 0.000000+0          2          2         20          09440 2151\n"
          " 4.633000+0 2.500000+0 1.000000+0 5.744900-4 3.100000-2 0.000000+09440 2151\n"
          " 1.000000-9 1.000000-9 1.000000-9 1.000000-9 1.000000-9 1.000000-99440 2151\n"
          " 1.000000-9 1.000000-9 1.000000-9 1.000000-9 1.000000-9 1.000000-99440 2151\n"
          " 1.000000-9 1.000000-9                                            9440 2151\n"
          " 0.000000+0 0.000000+0          0          0          0          09440 2  0\n"
          " 0.000000+0 0.000000+0          0          0          0          09440 0  0\n";
        auto begin = ENDF.begin();
        auto end = ENDF.end();
        long lineNumber = 0;

        int MAT = 9440;
        int MF = 2;
        int MT = 151;

        THEN( "the parameters can be extracted and verified" ){
          Isotope iso( begin, end, lineNumber, MAT, MF, MT );
          CHECK_NOTHROW(
            std::get<
              UnresolvedEnergyDependentFissionWidths >
              ( iso.resonanceRanges().front().parameters() )
          );
        }
        WHEN( "LFW is invalid" ){
          njoy::Log::info( "Testing invalid LFW" );
          int MAT = 125;
          int MF = 2;
          int MT = 151;
          THEN( "an exception is thrown" ){
            std::string ENDF =
              " 1.001000+3 1.000000+0          0         -1          1          0 125 2151    2\n"
              " 1.000000-5 1.000000+5          2          1          0          0 125 2151    3\n";
            auto begin = ENDF.begin();
            auto end = ENDF.end();
            long lineNumber = 0;

            CHECK_THROWS(
              Isotope( begin, end, lineNumber, MAT, MF, MT ) );
          }

          THEN( "an exception is thrown" ){
            std::string ENDF =
              " 1.001000+3 1.000000+0          0          2          1          0 125 2151    2\n"
              " 1.000000-5 1.000000+5          2          1          0          0 125 2151    3\n";
            auto begin = ENDF.begin();
            auto end = ENDF.end();
            long lineNumber = 0;

            CHECK_THROWS(
              Isotope( begin, end, lineNumber, MAT, MF, MT ) );
          }
        }
      }
    }
    WHEN( "LRF=2, all energy-dependent parameters" ){
      std::string ENDF =
        // " 3.809000+4 8.913540+1          0          0          1          03843 2151    1\n"
        " 3.809000+4 1.000000+0          0          1          1          03843 2151    2\n"
        " 6.000000+3 1.000000+5          2          2          0          03843 2151    3\n"
        " 0.000000+0 6.795900-1          0          0          3          03843 2151    4\n"
        " 8.913540+1 0.000000+0          0          0          1          03843 2151    5\n"
        " 5.000000-1 0.000000+0          5          0         84         133843 2151    6\n"
        " 0.000000+0 0.000000+0 0.000000+0 1.000000+0 0.000000+0 0.000000+03843 2151    7\n"
        " 6.000000+3 1.240600+4 0.000000+0 4.590100-1 2.050000-1 0.000000+03843 2151    8\n"
        " 7.000000+3 1.240600+4 0.000000+0 4.590100-1 2.050000-1 0.000000+03843 2151    9\n"
        " 8.000000+3 1.240600+4 0.000000+0 4.590100-1 2.050000-1 0.000000+03843 2151   10\n"
        " 1.000000+4 1.231300+4 0.000000+0 4.555800-1 2.050000-1 0.000000+03843 2151   11\n"
        " 1.500000+4 1.226300+4 0.000000+0 4.537400-1 2.050000-1 0.000000+03843 2151   12\n"
        " 2.500000+4 1.212200+4 0.000000+0 4.485000-1 2.050000-1 0.000000+03843 2151   13\n"
        " 3.000000+4 1.205400+4 0.000000+0 4.460100-1 2.050000-1 0.000000+03843 2151   14\n"
        " 4.000000+4 1.191800+4 0.000000+0 4.409700-1 2.050000-1 0.000000+03843 2151   15\n"
        " 5.000000+4 1.181000+4 0.000000+0 4.369800-1 2.050000-1 0.000000+03843 2151   16\n"
        " 6.000000+4 1.113000+4 0.000000+0 4.118100-1 2.050000-1 0.000000+03843 2151   17\n"
        " 7.000000+4 1.077400+4 0.000000+0 3.986200-1 2.050000-1 0.000000+03843 2151   18\n"
        " 8.000000+4 1.039500+4 0.000000+0 3.846100-1 2.050000-1 0.000000+03843 2151   19\n"
        " 1.000000+5 9.669700+3 0.000000+0 3.577800-1 2.050000-1 0.000000+03843 2151   20\n"
        " 8.913540+1 0.000000+0          1          0          2          03843 2151   21\n"
        " 5.000000-1 0.000000+0          5          0         84         133843 2151   22\n"
        " 0.000000+0 0.000000+0 0.000000+0 1.000000+0 0.000000+0 0.000000+03843 2151   23\n"
        " 6.000000+3 1.240600+4 0.000000+0 6.723800+0 2.050000-1 0.000000+03843 2151   24\n"
        " 7.000000+3 1.240600+4 0.000000+0 6.723800+0 2.050000-1 0.000000+03843 2151   25\n"
        " 8.000000+3 1.240600+4 0.000000+0 6.723800+0 2.050000-1 0.000000+03843 2151   26\n"
        " 1.000000+4 1.231300+4 0.000000+0 6.673600+0 2.050000-1 0.000000+03843 2151   27\n"
        " 1.500000+4 1.226300+4 0.000000+0 6.646600+0 2.050000-1 0.000000+03843 2151   28\n"
        " 2.500000+4 1.212200+4 0.000000+0 6.569900+0 2.050000-1 0.000000+03843 2151   29\n"
        " 3.000000+4 1.205400+4 0.000000+0 6.533500+0 2.050000-1 0.000000+03843 2151   30\n"
        " 4.000000+4 1.191800+4 0.000000+0 6.459600+0 2.050000-1 0.000000+03843 2151   31\n"
        " 5.000000+4 1.181000+4 0.000000+0 6.401200+0 2.050000-1 0.000000+03843 2151   32\n"
        " 6.000000+4 1.113000+4 0.000000+0 6.032400+0 2.050000-1 0.000000+03843 2151   33\n"
        " 7.000000+4 1.077400+4 0.000000+0 5.839300+0 2.050000-1 0.000000+03843 2151   34\n"
        " 8.000000+4 1.039500+4 0.000000+0 5.634100+0 2.050000-1 0.000000+03843 2151   35\n"
        " 1.000000+5 9.669700+3 0.000000+0 5.241000+0 2.050000-1 0.000000+03843 2151   36\n"
        " 1.500000+0 0.000000+0          5          0         84         133843 2151   37\n"
        " 0.000000+0 0.000000+0 0.000000+0 1.000000+0 0.000000+0 0.000000+03843 2151   38\n"
        " 6.000000+3 6.202800+3 0.000000+0 3.361900+0 2.050000-1 0.000000+03843 2151   39\n"
        " 7.000000+3 6.202800+3 0.000000+0 3.361900+0 2.050000-1 0.000000+03843 2151   40\n"
        " 8.000000+3 6.202800+3 0.000000+0 3.361900+0 2.050000-1 0.000000+03843 2151   41\n"
        " 1.000000+4 6.156500+3 0.000000+0 3.336800+0 2.050000-1 0.000000+03843 2151   42\n"
        " 1.500000+4 6.131600+3 0.000000+0 3.323300+0 2.050000-1 0.000000+03843 2151   43\n"
        " 2.500000+4 6.060800+3 0.000000+0 3.285000+0 2.050000-1 0.000000+03843 2151   44\n"
        " 3.000000+4 6.027200+3 0.000000+0 3.266800+0 2.050000-1 0.000000+03843 2151   45\n"
        " 4.000000+4 5.959100+3 0.000000+0 3.229800+0 2.050000-1 0.000000+03843 2151   46\n"
        " 5.000000+4 5.905100+3 0.000000+0 3.200600+0 2.050000-1 0.000000+03843 2151   47\n"
        " 6.000000+4 5.564900+3 0.000000+0 3.016200+0 2.050000-1 0.000000+03843 2151   48\n"
        " 7.000000+4 5.386800+3 0.000000+0 2.919600+0 2.050000-1 0.000000+03843 2151   49\n"
        " 8.000000+4 5.197500+3 0.000000+0 2.817000+0 2.050000-1 0.000000+03843 2151   50\n"
        " 1.000000+5 4.834900+3 0.000000+0 2.620500+0 2.050000-1 0.000000+03843 2151   51\n"
        " 8.913540+1 0.000000+0          2          0          2          03843 2151   52\n"
        " 1.500000+0 0.000000+0          5          0         84         133843 2151   53\n"
        " 0.000000+0 0.000000+0 0.000000+0 1.000000+0 0.000000+0 0.000000+03843 2151   54\n"
        " 6.000000+3 6.202800+3 0.000000+0 2.233000-1 2.050000-1 0.000000+03843 2151   55\n"
        " 7.000000+3 6.202800+3 0.000000+0 2.233000-1 2.050000-1 0.000000+03843 2151   56\n"
        " 8.000000+3 6.202800+3 0.000000+0 2.233000-1 2.050000-1 0.000000+03843 2151   57\n"
        " 1.000000+4 6.156500+3 0.000000+0 2.216300-1 2.050000-1 0.000000+03843 2151   58\n"
        " 1.500000+4 6.131600+3 0.000000+0 2.207400-1 2.050000-1 0.000000+03843 2151   59\n"
        " 2.500000+4 6.060800+3 0.000000+0 2.181900-1 2.050000-1 0.000000+03843 2151   60\n"
        " 3.000000+4 6.027200+3 0.000000+0 2.169800-1 2.050000-1 0.000000+03843 2151   61\n"
        " 4.000000+4 5.959100+3 0.000000+0 2.145300-1 2.050000-1 0.000000+03843 2151   62\n"
        " 5.000000+4 5.905100+3 0.000000+0 2.125900-1 2.050000-1 0.000000+03843 2151   63\n"
        " 6.000000+4 5.564900+3 0.000000+0 2.003400-1 2.050000-1 0.000000+03843 2151   64\n"
        " 7.000000+4 5.386800+3 0.000000+0 1.939200-1 2.050000-1 0.000000+03843 2151   65\n"
        " 8.000000+4 5.197500+3 0.000000+0 1.871100-1 2.050000-1 0.000000+03843 2151   66\n"
        " 1.000000+5 4.834900+3 0.000000+0 1.740500-1 2.050000-1 0.000000+03843 2151   67\n"
        " 2.500000+0 0.000000+0          5          0         84         133843 2151   68\n"
        " 0.000000+0 0.000000+0 0.000000+0 1.000000+0 0.000000+0 0.000000+03843 2151   69\n"
        " 6.000000+3 4.135200+3 0.000000+0 1.488700-1 2.050000-1 0.000000+03843 2151   70\n"
        " 7.000000+3 4.135200+3 0.000000+0 1.488700-1 2.050000-1 0.000000+03843 2151   71\n"
        " 8.000000+3 4.135200+3 0.000000+0 1.488700-1 2.050000-1 0.000000+03843 2151   72\n"
        " 1.000000+4 4.104300+3 0.000000+0 1.477600-1 2.050000-1 0.000000+03843 2151   73\n"
        " 1.500000+4 4.087700+3 0.000000+0 1.471600-1 2.050000-1 0.000000+03843 2151   74\n"
        " 2.500000+4 4.040500+3 0.000000+0 1.454600-1 2.050000-1 0.000000+03843 2151   75\n"
        " 3.000000+4 4.018100+3 0.000000+0 1.446500-1 2.050000-1 0.000000+03843 2151   76\n"
        " 4.000000+4 3.972700+3 0.000000+0 1.430200-1 2.050000-1 0.000000+03843 2151   77\n"
        " 5.000000+4 3.936800+3 0.000000+0 1.417200-1 2.050000-1 0.000000+03843 2151   78\n"
        " 6.000000+4 3.710000+3 0.000000+0 1.335600-1 2.050000-1 0.000000+03843 2151   79\n"
        " 7.000000+4 3.591200+3 0.000000+0 1.292800-1 2.050000-1 0.000000+03843 2151   80\n"
        " 8.000000+4 3.465000+3 0.000000+0 1.247400-1 2.050000-1 0.000000+03843 2151   81\n"
        " 1.000000+5 3.223200+3 0.000000+0 1.160400-1 2.050000-1 0.000000+03843 2151   82\n"
        " 0.000000+0 0.000000+0          0          0          0          03843 2  099999\n"
        " 0.000000+0 0.000000+0          0          0          0          03843 0  0    0\n";
      auto begin = ENDF.begin();
      auto end = ENDF.end();
      long lineNumber = 0;

      int MAT = 3843;
      int MF = 2;
      int MT = 151;

      THEN( "no exception is thrown" ){
        Isotope iso( begin, end, lineNumber, MAT, MF, MT );
        CHECK( 1 == iso.resonanceRanges().size() );
        CHECK_NOTHROW( std::get<
                         UnresolvedEnergyDependent >
                         ( iso.resonanceRanges().front().parameters() ) );
      }
    }
    WHEN( "LRF is invalid" ){
      int MAT = 125;
      int MF = 2;
      int MT = 151;
      THEN( "an exception is thrown" ){
        std::string ENDF =
          " 1.001000+3 1.000000+0          0          0          1          0 125 2151    2\n"
          " 1.000000-5 1.000000+5          2         -1          0          0 125 2151    3\n";
        auto begin = ENDF.begin();
        auto end = ENDF.end();
        long lineNumber = 0;

        CHECK_THROWS(
          Isotope( begin, end, lineNumber, MAT, MF, MT )
        );
      }

      THEN( "an exception is thrown" ){
        std::string ENDF =
          " 1.001000+3 1.000000+0          0          0          1          0 125 2151    2\n"
          " 1.000000-5 1.000000+5          2          3          0          0 125 2151    3\n";
        auto begin = ENDF.begin();
        auto end = ENDF.end();
        long lineNumber = 0;

        CHECK_THROWS(
          Isotope( begin, end, lineNumber, MAT, MF, MT ) );
      }
    }
  }
}
